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Abstract 

We study the processes of particle acceleration which take place in 
the field of a rotating black hole as part of a mechanism of formation 
of galactic jets within the first parsec from the central source, where 
gravitation is supposed to be dominant. We find the Lorentz factor 
that a stream of particles acquires as function of distance, when the 
orbital parameters vary slightly due to a local electromagnetic field or 
a pressure gradient. 

1 Introduction 

Jets emerging from active galactic nuclei are highly collimated structures (as 
revealed by radio maps) probably made of electron-positron plasma, which 
propagate in the intergalactic medium with relativistic velocities along most 
of their length (as implied by the detection of superluminal motion). It 
is then clear that one has to search for mechanisms which allow for both 
collimation and acceleration. 
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On large scales (typically from parsec to kiloparsec) jets are currently 
studied within the framework of magnetohydrodynamics, to which we owe 
most of our knowledge about these structures. Through analytical and nu- 
merical models, we are now able to say something about electrodynamic 
confinement of axisymmetric flows, electromagnetic extraction of energy 
from rotating black holes through force-free magnetospheres and asymptotic 
poloidal velocities of hydromagnetic flows, to mention but a few. What is 
shared by most of these models is the need for plasma injection into a rapidly 
rotating magnetosphere from below. 

Attention is then shifted to the subparsec scale, just outside the event 
horizon, where a primary acceleration mechanism must be at work. Recent 
developments (Sikora et al. 1996, O) have shown that such a mechanism 
cannot be powered by the radiation pressure of the accretion disk; this pres- 
sure, instead, causes the flow in the jet to decelerate by virtue of inverse 
Compton scattering with the plasma of the jet itself, with a maximum ef- 
ficiency when the plasma's Lorentz factor has reached values higher than 
an equilibrium ^ eq ~ 4. That is to say, jets, accelerated in the subpar- 
sec region up to 7j > 5, as revealed by VLBI measurements of superluminal 
motion in extragalactic radio sources, cannot avoid radiation drag. Notwith- 
standing this, we expect that at subparsec scale the gravitational field of a 
10 8 4- 1O U M black hole should still play a major role in determining par- 
ticles' motion; indeed the behaviour of individual particles is also that of 
the bulk of fluid elements in the guiding centre approximation. We then 
studied the combined effects of gravity and external physical perturbations. 
In de Felice and Carlotto (1997, hereafter Paper I) the collimating be- 
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haviour of geodesic orbits in the presence of a constrained variation of their 
energy and angular momentum was considered. The constraints are those 
which allow a particle, initially on a geodesic, to move on a nearby geodesic 
characterized by slightly varied parameters. In Paper I, this requirement was 
termed geodesicity condition. Following this line of thought, Karas and Dov- 
ciak (1997, |J) have estimated the rate of change of the orbital parameters 
of individual particles and have integrated this rate over a power law distri- 
bution of particles' energy. Their results confirm our claim in Paper I that 
the approximation of geodesic motion in the presence of small perturbations, 
is appropriate for modelling the primary collimation of a jet. 

Here we shall investigate whether the geodesicity conditions considered 
in Paper I, are compatible with the local Lorentz factor which is observed in 
galactic jets. Indeed we show that, provided there is a fine tuning between 
the stiffness of geodesic orbits and the effects of an external field, a large 
family of particle trajectories described by the Kerr metric not only collimate 
towards the axis of symmetry (see Paper I), but also accelerate, reaching 
values of the Lorentz factor 7, as measured by a local static observer, which 
are consistent with observations. Typically, we find 7 < 10 at lpc from the 
centre. 

In Section 2 we summarize the general relativistic collimation effect dis- 
cussed in Paper I, then in Section 3 we analyse the acceleration which test 
particles acquire under the condition of the mentioned collimation process. 
Behaviours of the local Lorentz factor 7 are found as functions of distance 
from the central source for both the cases of Lorentz forces arising from a 
local electromagnetic field and from pressure gradients. Comparison with 
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observations must be handled with care: some comments on this problem 
are made in Section 4. Finally in Section 5, we pay particular attention to 
the limits imposed by the geodesicity conditions on the distance scale where 
the collimation mechanism is allowed to operate. 

2 Geometrically induced collimation 

In Paper I, de Felice and Carlotto studied the tendency of vortical geodesies, 
in a Kerr background geometry, to collimate towards the axis of symmetry 
under a constrained variation of the constants of the motion. Their claim 
was that such a property might be astrophysically relevant to allow for jets 
primary collimation very close to the central black hole. Here we shall 
summarize that reasoning. 

Let the space-time be described by the Kerr metric in Boyer and Lindquist 
coordinates (t,r,9,ip): 

2Mr 1A 
ds 2 = -(1 ^)dt 2 - —u> sin 2 6 dtdtp + 

+-|dr 2 + £dfl 2 + ^sin 2 # dip 2 (1) 

where 

£ = r 2 + a 2 cos 2 # (2) 

A = r 2 + a 2 - 2Mr (3) 

A = (r 2 + a 2 ) 2 -a 2 A sin 2 9 (4) 
2Mar 

U = — (5) 
M and a being respectively the total mass energy of the metric source and 
its specific angular momentum (a = J/cM), both expressed in geometrized 



units (c = G = 1). 

Attention is focused on a particular family of geodesies, namely the vor- 
tical ones. These curves are gravitationally unbound (open orbits) and are 
characterized by the following conditions Q 

T>0 -a 2 T<L<a 2 T L < I 2 < {L ± fj? (6) 

4a 2 r 

In the absence of any external perturbation, I and E = y/T + 1 are constants 
of the motion and express, respectively, the azimuthal angular momentum 
(in units of uc) and the total energy (in units of uc 2 ) of the particle along 
the orbit. L is the separation constant of the Hamilton- Jacobi equation in 
the Kerr metric and is related to the square of the total angular momentum 
of the particle (de Felice, 1980, j§; de Felice and Preti 1998, §). 

The permitted values of 9 for the geodesic motion are confined (see figure 
[lj) to the area below the functional curves: 

l 2 (9,L,F) = sin 2 9(L + Ta 2 cos 2 9) (7) 

in the (I 2 , 9) space (for any fixed pair of L and T), along which 9 = 0. Given 
a constant value of I 2 in accordance with condition (||), vortical motion is 
found to be latitudinally confined within the range [0±, #2] determined by the 
intersection of the straight line I 2 = const with the functional curves (|7|). 
That is the reason why, vortical geodesies, which never cross the equatorial 
plane, are the most likely to leave the innermost part of an accretion disk 
surrounding a rotating black hole, through a spiralling motion. 

Collimation is studied with respect to the opening angle of a particle 

beam centred on the axis (9 = 0), taken for convenience to be equal to the 

1 Vortical geodesies are also those with r = I = L — 0, but we shall not consider them 
here. 
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angle which makes I 2 vanish: 

cos 2 ^ = -^. (8) 

The cornerstone of the reasoning is to presume that, in the innermost 
region of the field, relation (§) continues to hold even under a slow variation 
of the parameters T and L induced by some sort of external perturbation to 
the plain geodesic motion. 

In order to prevent the geodesic character of the motion being lost be- 
cause of the perturbations, we will impose "geodesicity conditions" with the 
effect of forcing each particle on a vortical geodesic to drift onto a nearby 
geodesic of the same type. 

In this case, differentiation of (||) with respect to the proper time r along 
the trajectories, leads to the basic equation of jet dynamics: 

„ d9n 9 o „ dT dL , . 

- 2Ltan6»o-r^ = a 2 cos 2 d — + — . 9) 
ar ar ar 

The existence of vortical trajectories in the space-time of a Kerr black 
hole, is a general relativistic effect which stems from the rotational properties 
of the metric, (see also O'Neil (1995) ||). While the vortical character of the 
orbits with parameters as in (6), is a natural consequence of gravitational 
dragging, the trend to axial collimation, as a result of a small perturbation, 
was quite unexpected. This effect is entirely due to the first term on the 
right-hand-side of equation (9) which contains the rotational parameter a. 
Evidently, when a = 0, there are no vortical orbits and no axial collimation; 
this implies that sufficiently small values of a would make the effect negli- 
gibly small. However, if we consider, as source of the orbital perturbations, 
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physical fields that share the same symmetries as the metric source, as ex- 
pected very near to the black hole, the collimation follows laws which do not 
depend explicitely on a (see equations (10) to (12) below), suggesting a kind 
of contradiction. Indeed, the dependence on a comes implicitly through the 
very form of the equations of motion and the equations of the perturbing 
fields. Parameter a fixes, from relations (6), the range of the permitted 
vortical orbits; when a becomes small, this range shrinks (see figure 1), so 
what decreases is not the amount of collimation, but the number of orbits 
which are involved. For this reason, this general relativistic effect may not 
be negligible at astrophysical scales. 

The explicit form of the variations (dT/dr) and (dL/dr) of the orbital 
parameters, depends on the nature of the perturbation we consider. As 
shown in Paper I, in the presence of a local electromagnetic field or of a 
pressure gradient, which are the two perturbations we are going to deal 
with, relations exist which link the coordinate to the energy parameter T. 
In case of energy gain (^7 > 0), such relations give rise to collimation laws, 
which we now recall: 

• Case of poloidal electromagnetic field (see section 3.1): 

— cospiralling 8 = const, orbits: 



sin a = sin I 



y/1 + i/r + 1 + l/(2r) 
V 1 + l/Fi + 1 + 1/(2I\) 

— counterspiralling 8 = const, orbits: 



1/4 



(10) 



sin 8 = sin ( 



r t (^t + r, + i/2 + r t ) 
r( v / T2Tr + i/2 + r) 



1/4 



(ii) 
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Case of pressure gradients (see section 3.2): 
— Either case of 9 = const, orbits: 

T^r + i)! 1 / 4 



sin V = sin ( 



(12) 



+ 1) 

where 6i and Tj are the initial values. In figure | and | we show the collimat- 
ing bahaviour of (|l0|), (|Tl| ) and (12) under the influence of the corresponding 



external perturbation which causes an increase of particle's energy. We have 
chosen as initial values, Tj = 0.1 at 0i = 7r/4. 

Distinction has been made between co-rotating (with the metric source) 
and counter-rotating orbits, collimation being stronger for the latter 0. Let 
us underline the physical meaning of these relations: geometry induced col- 
limation, in the presence of an electromagnetic field or a pressure gradient, 
occurs mainly when particles, on initially vortical geodesies, increase their 
energy with respect to infinity. However this is not the same as saying that 
they are locally accelerating, as we are going to see. 

3 Test particle acceleration 

Relations ( JToj ) - (11) - (p^) pressupose a knowledge of how V varies along 
the perturbed geodesic under the two kinds of external perturbations we 
are here considering. In this way, besides describing the collimation of the 
vortical geodesies, we will also deduce the behaviour of the local Lorentz 
factor 7 of the particles which leave the neighbourhood of the rotating black 
hole. 



2 Curves of figure 2 correct figure 4 of Paper I where they have been erroneously crossed. 
In that figure the two curves had to be considered independently and not to be compared. 
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Let us introduce a field of local static observers defined by a four- velocity: 

u=(-g 00 r 1/2 dt = e^d t (13) 

and let k be the particle four-velocity. The relation between the total par- 
ticle energy E, as measured at asymptotic distances, and the fundamental 
astrophysical quantity 7 locally measured by the observer u, is then given 
by: 

7 = -(u\k) = -e*(dt\k) =e^E (14) 
that is: 

7 



\ 



1 



r + i 



~a* COS^ tl 

where M, as stated, is the mass of the black hole. 

For the sake of comparison, let us first consider the behaviour of 7 for 
a particle in strictly geodesic motion, namely in the absence of external 
perturbations. 

Suppose the particle is moving outwardly on a geodesic with 6 = const = 
0.1 rad, arriving at infinity with 7oo = e^(°°^E = E = 1.1, starting from 
fj = rj/M = 1.5 with 7j = 1.1 e^' 6 ^ = 4.07 0. As we can see from fig[|, the 
particle is progressively slowed down as seen by a local static observer as it 
moves on its outwardly path. Such a trend reflects the attractive character 
of gravity, according to intuition. We have chosen a = a/M = 0.9981. 

In what follows we shall analyse the behaviour of 7 in a stream of parti- 
cles under the influence of an electromagnetic field and a pressure gradient, 
constrained however by the geodesicity conditions which guarantee the si- 
multaneous occurrence of collimation. 



3 Being IMq = 2.2 ■ 10 3 m in geometrized units, rt — 1.5 corresponds to a distance 
from a 10 8 M© black hole of n = n ■ M = 3.3 • lO^m ~ 10~ 5 pc. 



3.1 Effects of the electromagnetic field 

The electromagnetic field we consider arises locally from the potential: 

Or 

A = --^(dt-asin^edip) (16) 

where Q is the total electric charge. 

A charged particle of rest mass fi and charge q will deviate from geodesic 
motion by the term (neglecting radiation reaction): 

k r V r h = -Fijk j (17) 
A* 

where = 2duAji and k is, as stated, the 4-velocity of the particle. Since 
ko = —E, we easly deduce: 



dE qQ 



dr /iS 2 



(a 2 cos z 9 -r 2 )— + ra 2 sm 29 



(18) 



dr dr_ 

which tells us how E varies along the perturbed geodesic. Since the variation 
of E is much more sensitive to the variation of the coordinate r than to 9, 
decreasing as (M/r) 2 in the first case and as (M/r) 4 in the latter, we can 
take as full variation of E the following: 

dE qQ , 2 2 2 



dr /j,T, 2 



(r 2 - a 2 cos 2 0) (19) 



Recalling that E = \/T + 1 and using normalized quantities, relation (19) 
can be written as: 

8T - r 2 - a 2 cos 2 6 *- . , 

-5= = 2C T^ =5 9^vTTT 20 

or [r z + a z cos z 9) z 

where we have put C = -jjjj. 

This is the partial differential equation we were looking for. It can be 
coupled to the laws of collimation (|i~C|) and (|Tl|), providing a system in 
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the unknowns T(f) and 9(T(f)). This allows us to reduce (f2C|) from the 
form dfT = f(T, f, a cos 6) to the form dfT = g(T, f), exploiting the laws of 
collimation we have already found. The new partial differential equation can 
be numerically solved for r(f), telling us how T varies along the perturbed 
geodesic as a function of distance. On the other hand, knowledge of T(f) 
allows us to find the value of 6 as a function of r. Getting the Lorentz 



factor 7 = 7(f) would be the last step, through relation (15) with the laws 
of collimation to be used again. 

All seems to be quite smooth, except for the presence of parameter C: in 
Paper I it was shown that having a magnetic field B ~ 10 _3 G at the distance 
of lpc from a 10 s M@ black hole along the axis, implies Q ~ 2.3 • 10 11 m in 
geometrized units, so that Q/M ~ 1. Consequently, C is mainly the ratio 
q/l-t, which depends critically on the ionization degree. 

Since we are interested in the bulk motion of the material, rather than 
the motion of individual particles, we can accept the guiding centre approx- 
imation and look at our particle as a small cloud of hydrogen, say, whose 
specific charge we need to estimate. This requires taking into account the 
effects of photoionization due to the radiation field of the central source, the 
degree of recombination, the optical depth within the cloudlet itself and the 
gravitational redshift of the ionizing radiation. But first of all, we ought 
to know what kind of conditions on ratio q/fi comes from the geodesicity 
condition, so essential to our discussion. As stated in Paper I, the geodesic 
character of the motion can be approximately saved if orbital parameters 
vary slowly in time, geodesicity being better preserved where gravity dom- 
inates. This led the authors to evaluate the changes of energy of an orbit 
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that deviates from a geodesic under the effects of an electromagnetic field, 
written in terms of coordinates proper values variations (relations (41) of 
Paper I) as: 

dE . 
or 

q Q . , r 2 — a 2 cos 2 9 



5E\ 



6=cost 



-M- 



-,5r 



fiM (r 2 + a 2 cos 2 9) 2 

q_Q_fM_\ 2 f^ 1 ' 2 l-(a 2 /r 2 )cos 2 < 
fiM \ r ) 



(21) 



.Sy [1 + (a 2 /r 2 ) cos 2 #] 2 K M' 
where 5l r = (grr) 1 ^ 2 ^ is the proper-radial length as it would be measured 
by the particle itself. Since a variation of l r goes on with dynamical time 
while a variation of the parameter E goes with the perturbation time, the 
geodesicity conditions require that 5E < 5l r . 

More precisely, what we ask is that the timescale of a significant vari- 
ation of the physical parameters, T var , be longer than the dynamical time 
associated with a geodesic trajectory, r^ yn . Just for an order of magnitude 
estimate, let us take, as a significant dynamical time, the proper time a 
particle takes to reach the r = disk from the outer horizon on a parabolic 
trajectory, as was shown in Paper I: 



T~dyn 

On the other hand: 



2M 

1T 



i 



1 - 1 1? 



2M 
"3" 



E 



(22) 



(23) 



drj ^(r 2 - a 2 cos 2 0)f 

If we now recover from (Al) of Paper I the r component of the four-vector 

k tangent to a 9 = geodesic in which we choose, for simplicity, T ~ 1: 

k r _ ( r 2 + a 2 + 2Mr\ 1/2 
I r 2 + a 2 J 



(24) 
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we get: 

+ a 2 ? 12 (v ~ 

T V ar ~ ^ Q(r2 _ + ^ + 2Mr) l/ 2 I °J 

Therefore the geodesicity condition translated into timescale language re- 
quires that: 

I <r 3 /o g ± » 2 ) 5/2 = r rwn 

^ < 2 V V-a 2 )(f2 + a2 + 2f )i/2 ~ 1 j 

This is an explicit, though crude, expression of the geodesicity condition 

in the presence of the particular elettromagnetic field we have considered. 

It can be either considered as a condition on the ratio q/fj, if we fix the 

distance r from the source, or as a condition on the distance scale if we 

can fix the ionization. Taking as mentioned Q/M ~ 1 and a/M ~ 0.9981, 

(p6|) tells us that q/fi must be of the order < 10 at short distances from 

the centre (r = 2M), while it can increase outwardly up to q/fj, < 10 10 at 

r ~ 10 5 M ~ lpc. 

The astrophysical implications of these requirements have been discussed 
in Paper I, therefore we shall here derive the behaviour of 7(7*) with the 
assumption that the conditions on qj \i are satisfied^. That is, we are now in 
a position to integrate equation (|20| ) and get the behaviour of 7(f) through 
relation (|l^) , provided we take C values in accordance with condition ( |26|) . 

To a first approach, we will fix C to some constant (and low) value, so as 
to reproduce the behaviour of an approximately neutral plasma. Secondly, 
we will consider C linearly rising from the centre towards the outer regions, 

so as to reproduce the behaviour of increasingly ionized plasma, as it could 

4 Since the order of magnitude of the ratio q/fj, for an electron is, in geometrized units, 
q/fj, ~ 10 21 , we can deduce that for the godesicity condition to be satisfied, we ought to 
deal with cloudlets which are almost neutral near the horizon, and increase their ionization 
as we move outwards. Evidently, we shall refer to q//j, as the average degree of ionization 
in a volume element, rather than to the specific charge carried by a single particle. 
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be under the effect of a ionizing flux of radiation. That such a trend satisfies 
condition ( pff ) can be appreciated in figure [|, where it is compared with the 
plot of C ms . 

Figures (6) to (9) contain the results of the first approach. They show the 
7(f) profiles obtained through the mechanism just outlined, where we have 
chosen C = 2.5, 5, 10 on a distance scale ranging from fj = 2.3 to f = 100, 
that is to say from almost outside the event horizon to r ~ 10 -3 pc from 
the centre. Each figure shows two curves, one for co-spiralling orbits (solid 
line), and the other for counter-spiralling ones (dashed line). 

Figures (10) and (11) contain the results of the second approach. They 
show the 7 profiles obtained with linearly rising C values, that is to say with 
C ~ (3f, with f3 = 0.8 and (3 = 1.1. The distance scale has been enlarged up 
to the first parsec from the central source. Only one curve has been drawn 
in this case, because the profiles for co-rotating and counter-rotating orbits 
superimpose at such distances. 

The following aspects can be underlined: 

• we are facing an acceleration mechanism, which increases both V and 
7 up to asymptotic adjustments in dependence of parameter C. We 
typically find 7 ~ 2 4- 4 at 10 _3 pc from the centre, and 7 ~ 10 at the 
distance of Ipc. For constant values of C, slightly higher Lorentz factors 
are reached on co-rotating orbits with respect to counter-rotating ones; 
the difference disappears on larger scales when C is supposed to rise 
linearly. 

• we can appreciate the distinction between parameter T, which is linked 
to energy E and which is always raising, and the Lorentz factor 7, 
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which is the quantity measured by the local observer. Contrary to T, 
7 shows a rapid decay within the first r ~ 10M from the centre (see 
figure |7|), where the effects of gravity are supposed to overwhelm those 
of the electromagnetic field. As a matter of fact, 7 profiles in this inner 
region can be almost superimposed on those we found for the geodesic 
motion, proving the stiffness of geodesic orbits very near the central 
black hole. 

• In our calculations we have chosen an initial value of Tj = 0.1 at 
fj = 2.3 and 9i = vr/4, meaning an initial value of ji ~ 2.3. What can 
power particles to such a relatively high value of the Lorentz factor just 
outside the event horizon? The question is still open: some version of 
the Penrose mechanism has been proposed (Reva K. 1995, ||]). 

• we should pay attention to the fact that the mechanism we have just 
explored coexists with radiation pressure effects. As we have men- 
tioned in the introduction, inverse Compton interaction of the plasma 
in the jet with the radiation field produces a net deceleration, which 
always arises when the Lorentz factor of the bulk motion exceeds an 
equilibrium regime given by j eq ~ 4, depending on the radiation field 
distribution one adopts. Though in this Paper we have not considered 
radiation field as a possible physical perturbation, we can argue, look- 
ing at our 7 profiles, that radiation drag will not occur within 10 -3 pc 
from the centre, where 7 does not exceed the value of 5. Far from the 
black hole, at lpc from it, we find 7 ~ 10, so there might be a per- 
ceivable effect of radiation drag tending to lower 7 to its equilibrium 
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value. But far from the black hole, according to Sikora et al. (1996), 
such an effect decreases in importance. 

3.2 Effects of the pressure gradient 

Let us now consider the effects on the geodesic motion of the perturbation 
represented by pressure gradients, which probably characterize the inner 
parts of an active galactic nucleus: 

P ~ P a (27) 
p ~ r~ n (28) 

with 1 < n < 3, as suggested by several polytropic models. 

These pressures and densities are assumed to be those of a perfect fluid, 
whose elements approximate the behaviour of the emerging particles. The 
relativistic Euler equation leads to a variation of the energy of each fluid 
element along the perturbed trajectory, given by: 

(IT p + pdT 

For the same reasons we appealed to in the electromagnetic field perturba- 
tion case, ( p9| ) transforms into: 

d r r = - 2(r + 1) d r p (30) 
p + p 



which is analogous to equation (|19|). It can be further developed if we take 
p = Ap a and p = Br~ n , without specific knowledge of the constants A and 

^Str + iy <■-->- 

1 + Pf n ( 1 - Q ) 
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"T A DQ-l 

where we nave put V = Mn(a _i) . 

Once again we need to evaluate energy changes on a 8 = const perturbed 
geodesic, in order to see what kind of geodesicity condition we have to 
account for. From Paper I we deduce: 



r dp 
p dr 



5E\e =cost < E 
with the inferred requirement that, at least: 

< l; 



M\ {A\ l / 2 x fl r 



(32) 



r dp 
p dr 



(33) 



such a condition is equivalent to the requirement that the sound speed in 
the comoving frame of the fluid is non-relativistic: 

1 



v s < 



n 



(34) 



The last step is now to choose realistic values for a and n. Theoretical 
arguments combined to semiempirical estimates of pressure based on VLBI 
measurements (see Begelman 1984, ||) seem to converge on a = 2, n = 2 
as the best choice^]. Therefore, the differential equation to integrate is, from 



dfT 



8£>(r + 1) 



r 3 + T>r 

with the associated geodesicity condition that, from (34), is: 

1 



v, < 



V2 



(35) 



(36) 



This latter condition implies a lower limit to the distance scale where the 

mechanism can operate. In order to see this, let us write down the behaviour 

5 a and n are not independent: for polytropic gases of the kind we are discussing we 
know that n — a /(a — 1). 
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of v s upon r as it results if we remember the definition vl = dp/ dp with our 
expressions for pressure and density: 

v 2 s ~ aW"^- 1 ) (37) 

and, with our choices for n and a: 

v s ~ v^f -1 (38) 

Now, given a value of P, (^) yields the behaviour of v s , which must satisfy 
(^). Hence there exists a lower limit in the distance scale, given by fib = 
2Vf>. 

We have numerically solved equation (|35|) using low values of T>. The 
Lorentz factor 7 comes from fll5|), with the appropriate law of collimation 



offered by (12). Figures (12), (13) and (14) show the results of our cal- 
culations: the distance scale extends from a point r, = 2vD (in all cases 
near the event horizon) to ~ 10~ 2 parsec from the centre. As we can see, 
the Lorentz factor shows a maximum at distances an order of magnitude 
farther than the event horizon from the centre, depending on parameter T>. 
At large distances, 7 tends to decrease, adjusting itself to a constant value. 
The mechanism seems to be most efficient in the very central regions of the 
field. 

4 Comparison with observations 

Some facts make it difficult to speak of an effective comparison with obser- 
vations: 

• what can be observed and measured, through the analysis of VLBI 
maps, is the proper motion of a radio pattern along the jet. Knowledge 
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of the redshift z of the source, combined with a hypothesis on the 
cosmological parameters qo and Hq, allow us to derive the apparent 
velocity f3 app at which the radio components are seen to be moving. 
Superluminal motion is detected over a wide distance range, starting 
from 10~ 2 pc from the core, up to 10 2 -f- 10 3 pc far away from it. 

The Lorentz factor corresponding to the superluminal motion that 
is usually observed with the above mentioned method is assumed to 
coincide with the Lorentz factor of the plasma bulk motion along the 
jet, responsible for the Doppler boosting of the radiation. 

• in such a manner, statistical analysis carried out on a sample of sources, 
for which superluminal motion has been detected, yield the following 
mean values of the Lorentz factor (taken from data of Ghisellini et al., 
1993, @) : 



Sources 


7 


BL Lac 


10.47 ±1.35 


CDQs 


16.98 ± 1.25 


LDQs 


14.45 ± 1.38 



Table 1: Mean values of the Lorentz factor for the sample of sources in the paper 
of Ghisellini et al. CDQs and LDQs refer respectively to core dominated quasars 
and lobe dominated quasars. 

What these data do not tell us is the distance from the source where 
the Lorentz factor has been measured, while our 7 values refer to dis- 
tances within the first parsec from the centre. 

• our 7 profiles have been found in the single-particle approach for the 



19 



case of the electromagnetic field perturbation, and in the hydrody- 
namical approach for the case of the pressure gradient perturbation. 
This guiding centre approximation is the first step towards a full mag- 
netohydrodynamical treatment. 

• Lorentz factors 7 derived from observations are found within the frame- 
work of special relativity, assuming relativistic motion of the outgoing 
plasma along a direction making an angle (p with respect to the line 
of sight. Therefore, we have stressed the importance of having intro- 
duced a local observer, because around him special relativity holds, so 
that what he measures is approximately what we measure from earth, 
apart from negligible problems due to the choice of a particular local 
observer (we found it useful to choose a local static observer). 

What is more, a direct comparison with observations is impossible unless we 
take into account the major role of large scale magnetic fields in extracting 
rotational energy from the black hole. Instead, what can be reasonably 
said is that the mechanisms we have proposed serve as primary acceleration 
mechanisms able to inject the plasma into the magnetosphere with values 
of the Lorentz factor 7 around 5 10. 

5 The distance scale 

The most important restrictions to the efficiency of the acceleration mech- 
anisms we have just outlined are represented by the geodesicity conditions, 
which turn out to estabilish conditions on the ionization degree, combined 
with conditions on the distance scale where the mechanisms can operate. 
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There does not seem to be any upper limit to the distance scale for the 
mechanisms we have developed, while there is a potential lower limit for both 
types of perturbation we have considered. Let us justify these statements 
starting from the electromagnetic field perturbation. 

If we come back to figure ||] we can infer the existence of a lower limit 
to the distance scale: since the curve describes the maximum degree of 
ionization allowed, a high constant value of q/(i would violate the geodesicity 
condition at short distances. The higher the value of <?///, the farther the 
lower limit is from the centre. That is the reason why in our calculations we 
have chosen low constant values for parameter C. Similarly, a linearly rising 
degree of ionization satisfies the geodesicity condition all along the distance 
scale, as the figure shows. 

On the other hand, it is the very nature of the electromagnetic potential 
to prevent the existence of an upper limit to the distance scale. In fact, 
both At and A v of potential ([!(]) vanish at infinity. This suggests that the 
effects of the electromagnetic field would vanish as we moved outwardly^] and 
we have found that the behaviour of V confirms this conclusion: V adjusts 
itself to a constant value far from the centre, as would be the case of a 
plain geodesic motion {5E — * 0). Only the case of a degree of ionization 
intersecting the limiting curve of figure @ from below would violate the 

geodesicity condition far from the centre, giving rise to an upper limit of the 

6 With a different choice of the electromagnetic field, such as that cited by Karas & 
Dovciak: 

A t = /3a[r£ _1 (f + cos 2 9) - 1] 

A v = /3sin 2 e[l/2(r 2 + a 2 ) -a 2 rE _1 (f + cos 2 0)] 

/3 ~ 5 • 1CT 8 (B /10 i Gauss)(M /fO 8 M©) ~ 10" 8 , the situation would have been different 
and we would have probably found an upper bound too. 
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distance scale, but that seems to be quite unphysical. 

Now let us turn to the pressure gradient case. We have already shown 
the existence of a lower limit in the distance scale, depending on the value 
of T>: the higher the value of T>, the farther the lower limit from the centre. 
No upper limit can be found either: that is physically suggested by the fact 
that pressure gradients vanish at infinity as dp/dr ~ r~ 5 so v s decreases 
monotonically in such a way that condition (|34] ) is always satisfied far from 
the black hole. 

In this sense, the acceleration mechanisms we have proposed need par- 
ticular care near the centre, where we have to hinder the effects of extreme 
physical conditions, while they die a natural death far from the black hole. 
Therefore, there is no point in extending the integration beyond the first 
parsec from the centre, since neither the local electromagnetic field, nor the 
pressure gradients have any more appreciable effect on the motion at those 
distances. 

6 Conclusions 

Radio components in jets emerging from active galactic nuclei are found 
to undergo superluminal expansion both near the core, at 10~ 2 pc from the 
centre, and far away from it, at 10 2 — 10 3 pc. Let us just cite 3C 273, that 
shows components in superluminal motion at ~ 50pc from the core (S.C. 
Unwin 1989, [fuj); 3C 345, studied at 22GHz, which shows superluminal 
expansion on the parsec scale (5 — 25pc from the core), and 0836+71, which 
shows superluminal motion at ~ 220pc from the centre. 

Therefore, observations suggest that jets are born with high Lorentz 



22 



factors. In Paper I we indicated that these structures emerge already colli- 
mated in the inner region of an AGN, instead of being collimated by large 
scale magnetic fields. In this Paper we have tried to complete the picture, 
testing the capability of a properly perturbed gravitational field to produce 
relatively large local Lorentz factors of escaping particles within the first 
parsec from the centre. 

We are now in a position to talk of collimation after acceleration, which 
would occur in the inner region of an AGN, where the effects of the black 
hole space time geometry mix with those of external perturbations, namely 
a local electromagnetic field and pressure gradients. Important restrictions 
on the ionization degree of a testing cloudlet have been introduced, in order 
to preserve the geodesic character of the motion; specifically, we have found 
that a degree of ionization which rises linearly with distance from the core 
is compatible with quasi-geodesicity. 

The nature of the interaction between the charged cloudlets and the 
electromagnetic field has been studied limitedly to the Lorentz forces which 
arise: more subtle effects such as Compton losses have been neglected. This 
idealized picture can be justified if we consider the order of magnitude of 
the competing forces, namely gravoinertial against Lorentz ones. As a con- 
sequence, the proposed mechanisms show most of their efficiency within the 
first parsec from the centre, producing values of the local Lorentz factor 7 
in the range 54- 10, depending on the parameters introduced. We claim that 
the nozzle appealed to by some hydrodynamical models could be a point in 
the region between 10 2 4- 10 3 r s from the core, where we find values of the 
Lorentz factor between 4 4 8 for the two cases of perturbation studied. 
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It is to the distance scale to which we refer when we define the above de- 
veloped acceleration mechanisms as primary acceleration mechanisms, and 
we can but confirm the need to look for different physical processes to ac- 
count for the velocity regime of jets, starting from the subparsec region and 
ending in the external radio lobes. 

Acknowledgments Thanks are due to Prof. Mary Evans Prosperi for 
correcting the English text. 
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Figure 2: Behaviour of the angle of collimation for corotating and counter- rotating 
outgoing particles which increase their energy with respect to infinity, under the 
influence of a poloidal magnetic field. 1^ = 0.1, 8i = n/A. 
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Figure 3: Behaviour of the angle of collimation for outgoing particles (corotating 
and counter rotating) which increase their energy with respect to infinity, under 
the influence of a pressure gradient. I\ = 0.1, Q% = 7r/4. 
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Figure 4: Lorentz factor for a particle on a plain geodesic motion. 9 = 
0.1 rod, E = 1.1 
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Figure 5: The figure shows how a straight line representing a linearly rising ion- 
ization (q/ fJ> — 4f) is always below the curve C ms in equation (26) which limits the 
validity of geodesicity condition. 
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Figure 7: Behaviour of the Lorentz factor 7 corresponding to the values of T 
represented in figure [| 
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Figure 8: Behaviour of the Lorentz factor 7 upon distance for co-rotating and 
counter-rotating orbits under the effect of the electromagnetic field. C — 5, Tj = 0.1, 
ft = 2.3. 
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Figure 9: Behaviour of the Lorentz factor 7 upon distance for co-rotating and 
counter-rotating orbits under the effect of the electromagnetic field. C = 10, I\ = 
0.1, n = 2.3. 
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Figure 10: Behaviour of the Lorentz factor 7 for co-rotating and counter-rotating 
orbits with a linearly rising degree of ionization: C = 0.8f; distance scale enlarged 
to the first parsec from the centre . 
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Figure 11: Behaviour of the Lorentz factor 7 for co-rotating and counter- rotating 
orbits with a linearly rising degree of ionization: C = l.lf; distance scale enlarged 
to the first parsec from the centre . 
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Figure 12: Behaviour of the Lorentz factor 7 for co-rotating and counter- rotating 
orbits under the effect of the pressure gradients. f> = 5, Tj = 0.1, 6i = 7r/4. 
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Figure 13: Behaviour of the Lorentz factor 7 for co-rotating and counter-rotating 
orbits under the effect of the pressure gradients. f> = 10, T, = 0.1, Oi = 7r/4. 
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Figure 14: Behaviour of the Lorentz factor 7 for co-rotating and counter-rotating 
orbits under the effect of the pressure gradients. f> = 15, T, = 0.1, Oi = 7r/4. 
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